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ABSTRACT 


Objective: To estimate the hydrological drought frequency for upper Tana River basin in Kenya using absolute Stream flow 
Drought Index (SDI) and modified Gumbel technique. The frequency of drought event of a defined severity for a defined return 
period is fundamental in planning, designing and operation of water storage systems in the basin. 

Materials and Methods: Based on a 41-year (1970-2010) stream flow data, hydrological droughts of 2, 5, 10, 20, 50, 100, 200, 
500 and 1000-year return periods are evaluated based on the stream flows, Stream flow Drought Index (SDI) and a simplified 
mathematical model for hydrological drought estimation which is formulated using Gumbel’s technique. 

Results: The absolute SDI increases while the magnitude of the stream flow decreases with return period. The minimum and 
maximum drought events were exhibited in gauge stations 4AC03 and 4CC03 with absolute SDI ranging from 0.667 to 1.265 and 
1.213 to 2.42, and corresponding stream flows of 4.341 to 2.719 and 18.246 to 1,021m 3 /s for a 2 and 1000-year return period 
respectively. 

Conclusion: A simplified mathematical model for estimating hydrological drought event that uses mean flows of the annual 
minimum and average of the first three minimum stream flows as input variables is formulated for different return periods for the 
river basin. 

Key Words: Upper tana River basin, SDI, Hydrological drought, Return period, Gumbel technique, Drought frequency, Math¬ 
ematical model 


INTRODUCTION 

Hydrological drought is a natural hazard associated with 
water deficiency in a hydrological system. Drought may 
be manifested in below average water availability such as 
stream flow in rivers, quantity of water in reservoirs, lakes 
and ground water (Tsakiris, 2009; Mishra and Singh, 2010; 
Sheffield and Wood 2011). Hydrological drought decreases 
the availability of water resources (Liu et al., 2012) in river 
basins, adversely impacting on economic aspects (Carrol et 
al., 2009; Van Vliet et al., 2012) social dimensions such as 
increased human conflicts and mortality rates (Garcia-Herre- 
ra et al., 2010) and ecological systems (Lake 2011, Lewis et 
al., 2011). There is need to understand the drought events in 
order to develop drought mitigation mechanisms in river ba¬ 


sins (Wambua et al, 2014). Hydrological drought impacts on 
large areas and large human population and may be triggered 
by climate change and /or variability (Mondal and Mujum- 
dar, 2015). Like other drought events, hydrological drought 
is considered to be a ‘creeping hazard’ because it develops 
slowly, it is not easily noticed, covers extensive areas and it 
lasts for long a period of time with adverse impact on eco¬ 
logical systems and socio-economic development (Liu et al., 
2015; Van-loon, 2015). In addition to hydrological droughts, 
other types of droughts Include meteorological agricultural/ 
soil moisture droughts and socio-economic drought. How¬ 
ever, according to Van-loon and Laaha (2015), hydrological 
drought has the most significant effects almost across differ¬ 
ent sectors as shown in Table 1. 
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Table 1: The major drought impacts for different drought categories 


Impact category 

Impact sub-category 

Hydrological drought 

Meteorological 

drought 

Agricultural/ 

soil moisture drought 

Agriculture 

Rain-fed 

X 

X 

X 


Irrigation 



X 

River basins/ecosystems 

Terrestrial 

X 

X 

X 


Cooling 

X 



Energy and industry 

Hydro-power 

X 




Cooling water 

X 



Navigation 


X 



Drinking water 


X 



Recreation 


X 




Source: Van-loon (2015) 


The key parameters of droughts are the longest duration and 
highest severity for a defined return period. Such parameters 
aid in designing water storage systems capable of withstand¬ 
ing effects of droughts (Kyambia and Mutua, 2014). Since 
occurrence of drought contributes to adverse socio-economic 
impacts, they need to be quantified so as to improvise coping 
and/or mitigation mechanise. Thus hydrological drought es¬ 
timation using stream flow data for a defined range of return 
period in a river basin is crucial. The commonly used return 
periods are 2, 5, 10, 20, 50, 100, 200, 500 and 1000 years. A 
50-year hydrological drought is defined as the drought mag¬ 
nitude which is equalled or exceeded, on average, once per 
50 years. 

The original version of Gumbel approach in the prediction 
requires computation of coefficient of variation (Cv) and 
determination of expected mean (y ) and standard deviation 
(rrj as given in the Gumbel’s table. However, in this paper, a 
mathematical model is formulated for hydrological drought 
estimation using annual minimum average stream flow and 
the mean of three lowest stream flows from the recorded 
data as the main input variables of a modified mathematical 
model. From an engineering point of view, design of water 
storage structures requires critical information of longest du¬ 
ration, largest severity for a specific return period. Extreme 
drought may be treated as a stochastic variable that is chal¬ 
lenging to estimate. For practical engineering work, extrapo¬ 
lation and interpolation of drought frequency is crucial for 
reliable designs. Gumbel (1958) put forth a method for esti¬ 
mating flood frequency. Such a technique may be modified 
and applied in drought frequency estimation. The Gumbel 
extreme value distribution for instance has been applied in 
drought studies in Greece (Dalezios et al., 2000) river basin 
and Dudhkumar River (Asad et al., 2013). 



Hydrological drought process 

The occurrence of hydrological droughts is considered to 
exhibit stochastic characteristics and thus complex in na¬ 
ture. Hydrological droughts are influenced by hydrological 
processes of the hydrologic cycle (Peters et al., 2006; Vidal 
et al., 2010) such as precipitation, evapotranspiration, soil- 
water storage, runoff flow on land and streams, and ground- 
water recharge or discharge. The fundamental cause of hy¬ 
drological drought is climatic change and/ or variability. For 
instance, an abnormally prolonged precipitation deficit leads 
to low input into hydrological system. Droughts may be trig¬ 
gered by anomalies in temperature in large scale atmospheric 
and or oceanic patterns and low sea temperature. For any 
river basin, the rate of depletion of soil-moisture is a func¬ 
tion of antecedent moisture condition, evaporation from bare 
soil surface, evapo-transpiration from vegetated areas, deep 
percolation of water into the groundwater and runoff into 
stream networks. For a dry season, runoff and drainage are 
significantly low while potential evapo-transpiration may be 
high as a result of increased solar radiation, vapour pressure 
deficit and wind velocity. During an extreme drought event, 
soil-moisture may be depleted to wilting point below which 
plants significantly undergo wilting due to response to the 
moisture decline in the soil media. This condition leads to 
reduction in actual evapo-transpiration and locally generated 
precipitation. The depletion of soil-moisture leads to de¬ 
crease in recharge of the groundwater storage. Soil-moisture 
is significantly influenced by the quantity of precipitation, 
recharge, discharge and aquifer storage and transmisivity 
characteristics. The relationship between precipitation, soil 
moisture, runoff, recharge, discharge, ground water and dis¬ 
charge is well explained using the hydrologic water balance 
Equation 1. 

R - P - ET - AS - G (1) 
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Where; 

R=Runoff (mm) 

^precipitation (mm) 

ET= E v apo-1 ran sp i rat ion (mm) 

/I.S-Ghange in soil-water storage (mm) 

G= Ground water (mm) 

The above relation presents an old concept of hydrology 
and has been well researched in numerous catchments in the 
world. However, the application of the hydrologic water bal¬ 
ance relationship in drought studies is a relatively new con¬ 
cept (Van loon, 2015). Any climate change and /or variabil¬ 
ity directly affect precipitation and evapo-transpiration, and 
indirectly influencing the runoff, soil-moisture storage, and 
groundwater components of the water balance model. Mete¬ 
orological droughts propagates into other types of droughts 
through processes of runoff, stream flow, recharge and dis¬ 
charge which are mainly influenced by river basin charac¬ 
teristics and climatic change or variability. However, the fre¬ 
quency of occurrence of drought is not well researched and 
thus not understood for numerous river basins in the world. 

Objective 

The objective of this research was to estimate hydrological 
drought frequency using absolute Stream flow Drought In¬ 


dex (SDI) and modified Gumbel’s technique for upper Tana 
River basin, Kenya. 


MATERIALS AND METHODS 

Description of upper Tana River basin 

The Tana River basin from which upper Tana River basin is 
delineated is the largest river basin in Kenya (Jacobs et al., 
2004; WRMA, 2010). It lies between latitudes 00° 05’ and 
01° 30’ south and longitudes 36° 20’ and 37° 60’ east. The 
upper Tana River basin has an area of 17,420 km 2 (Figure 1). 
The basin plays a critical role in regulating the hydrology of 
the entire basin (IFAD, 2012), and in the process, it controls 
the hydro-electric power generation within the Seven-Folk 
dams downstream of the Tana River. The basin is very criti¬ 
cal in Kenya as it drives the socio-economic development 
through water supply and agricultural production. The el¬ 
evation of the upper Tana River basin ranges from approxi¬ 
mately 730 m to 4,700 m above mean sea level (a.m.s.l.). 
These elevations are adjacent to Kindaruma hydro-power 
dam and Mount Kenya respectively. The dominant soil types 
in the basin are Andosols, Nitosols, Ferrasols and Vertisols 
at higher, middle and lower elevations respectively (Jacobs 
et al., 2004). 



Figure 1: The location of the upper Tana River basin in Kenya 
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Precipitation and temperature vary spatially across the entire 
river basin. The annual precipitation at Mount Kenya and the 
Aberdares ranges is 1800 mm (Otieno and Maingi, 2000). 
In the mid elevation of 1200 to 1800 m a.m.s.l., the annual 
rainfall ranges from 1000 to 1800 mm, while the lower el¬ 
evations at 1000 m, and receive annual rainfall of 700 mm. 
The basin is characterized by seasonal rainfall fluctuations 
as influenced by orographic forces (Saenyi, 2002). Subse¬ 
quently, this leads to seasonal variation of stream flows in 
Tana River. Generally the basin experiences bimodal rainfall 
pattern which is triggered by inter-tropical convergence zone 
(Wilschut, 2010). The two main rain seasons as shown in 
Figure 2, are distributed in the months of March to June, and 
September to December where the monthly average precipi¬ 
tation is considerably high compared to the other months. 



o e 

O cz 


Q= the probable hydrological drought discharge with a re¬ 
turn period of T years 

C =coefficient of variation 

v 

Q =The mean hydrological drought discharge (m 3 /s) 

K= frequency factor 

The frequency factor and the coefficient of variation are de¬ 
termined from the relation: 


K = 


(y T ~y„) 


(3) 



( 4 ) 


y T 


= - In In 


r t > 

J 1 T, 


(5) 


Where; 

v n , =expected mean 

o n ^standard deviation of reduced drought ex¬ 
tremes estimated from Gumbel’s Table 


In an attempt to simplify Equation (5), Powell (Asad et al., 
2013) developed an equation to estimate K using the relation: 


Figure 2: The monthly average precipitation 


K = 


(Vh) x (0.5772 + lnlnT / T - 1) 
n 


( 6 ) 


Data acquisition 

Stream flow data used in the present study was obtained 
from the Ministry of Water and Irrigation, and Water Re¬ 
sources Management Authority (WRMA) for eight stations 
for a period of 41 years (1970-2010). Data from eight sta¬ 
tions with consistent data that had less than 20% missing 
data was selected for the study. The Double mass curve was 
used to check for the data consistence. 

Gumbel’s extreme value (EV1) method 

Gumbel’s method was originally developed for flood esti¬ 
mation. However, it has previously been adopted in drought 
studies (Dalezios et al., 2000). The form of Gumbel tech¬ 
nique used for the present study for estimating extreme 
drought event is expressed as: 


Although the above limction improved the method of esti¬ 
mating the frequency factor K, which can now be computed 
based on return period T and not number of years of record, 
the method for calculating Cv still remains as suggested by 
Gumbel. 

Stream flow drought index 

A drought index is an integration of either one or more of 
hydro-meteorological variables such as precipitation, stream 
flow, soil moisture, temperature, ground water, water reser¬ 
voir volume or level (Sun et al, 2011). In this study. Stream 
flow drought index (SDI) that uses stream flow data is 
applied. The SDI for each gauged station was deter¬ 
mined using the following relation: 


Q T = Q{\ + KC,) 


Where; 


( 2) 

C7, 

Where; 


(7) 
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,S7.)/-strcam flow drought index for i th hydrological month 
Q-stream flow for the i tb hydrological month 
A=length of period of data record/reference period 


cr =the standard deviation of the cumulative stream 

k 

flow volumes for k lh reference period 

Table 2: Definition of states of drought based on SDI 


State 

Drought description 

Criterion 

0 

Non drought 

SDI > 0.0 

1 

Mild drought 

—1.0< SDI < 0.0 

2 

Moderate drought 

—1.5 < SDI < —1.0 

3 

Severity drought 

SDI < -2.0 

4 

Extreme drought 

SDI < -2.0 


The original function was developed by Gumbel (Gumbel, 
1958) for extreme flood estimation that used data that exhibit 
positive values. In this research, the stream flow drought in¬ 
dex with negative values as shown in Table 2 represent the 
period of drought episodes. These negative values are con¬ 
verted to their corresponding absolute values and fitted to 
Equation (2). Then the SD/for 2, 5,10, 20, 50, 100,200, 500 
and 1000-year return periods were computed using Equa¬ 
tion (7). The resulting SDI data was arranged into ascending 
order alongside the corresponding stream flow. The stream 
flow corresponding to the computed SDI of rank m for spe¬ 
cific return period was selected. Those stream flow values 
without corresponding SDI were interpolated using the rela¬ 
tions: 

Q m = Qo + soi Zf DI x(SDI '" ~ SDI ^ < 8) 

Where; 

Q m = the stream flow of rank m and specific return period 
(m 3 /s) 

gj=higher rank stream flow (m 3 /s) 

Q Hower rank stream flow (m 3 /s) 

SD/=the higher stream flow drought index 
SDI= the lower rank stream flow drought index 
SDI m = interpolated value of stream flow drought index 


Table 3. Absolute stream flow drought index (SDI) 
for computation of Q m 


Absolute SDI 

Stream flow Q m (m 3 /s) 

0.4609 

4.8339 

0.4809 

4.7814 

0.5229 

4.6708 

0.5692 

4.5494 

0.5753 

4.5334 

0.5797 

4.5217 

0.5873 

4.5018 

0.6119 

4.4371 

0.6485 

4.3408 

0.6790 

4.2610 

0.7226 

4.1463 

0.7246 

4.14098 


RESULTS AND DISCUSSION 

For instance for a 2-year return period, whose computed 
SDI is 0.6665, Q value was interpolated using Equation 
(8) based on the data in Table 3. 

The results of the fitted curves show that the absolute SDI in¬ 
creases with the return period in all gauged stations (4AB05, 
4BC02, 4AC03 and 4AD01) as given in Figure 3. For in¬ 
stance, hydrological droughts represented by magnitude of 
absolute SDI of 0.667and 1.265 are equaled or exceeded 
once on average every 2 and 1000 years respectively. The 
same applies to the other hydrological droughts of defined 
absolute SDI. For water resources managers, data on stream 
flow is important for ease of water resources planning and 
management. Thus in this research, the absolute SDIs are 
tied to their respective stream flow magnitudes. 



-•-SDI -*-Q 

a) 
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Using Table 2 and results from Figure 3, the critical points 
are identified. The critical point is the level of hydrological 
drought beyond which the water facilities is significantly af¬ 
fected by drought. For instance the critical point for 4AB05 
(Amboni) gauge station, as shown by dotted line, coincides 
with return period of 28 years with absolute SDI of 0.92 and 
stream flow of 3.6 m 3 /s (Figure 3a), while that of 4BC02 
(Tana sagana) is 20 years with absolute SDI of 1.2 and stream 
flow of 45 m 3 /s (Figure 3b). In this case, if a water resource 
system is to be designed for example at the gauging stations 
4AB05 and 4BC02, the systems should be designed for re¬ 
turn periods of less or equal to 28 and 20 years respectively. 

Form Figure 4, the results show that the ratio of Q T /Q rep¬ 
resented by Y increases with C for different return periods. 
This confirms that the Gumbel method is also applicable in 
drought frequency estimation just like in flood frequency 
analysis (Al-Mashindani etal., 1978). 


<0 



—•—SDI Q 

d) 

Figure 3: The relationship between the Q m , SDI and return 
period for (a) 4AB05 (Amboni), b) 4BC02 (Tana Sagana) c) 
4AC03 (Sagana) and (d) 4AD01 (Gura) gauge stations 


Development of the modified mathematical 
model 

Gumbel method has been modified before for flood studies. 
However, scanty research as far as its application in drought 
studies is concerned. In this research the principles used by 
Al-Mashindani (1978) in flood assessment was used in hy¬ 
drological drought estimation for upper Tana River basin. 
The value of Q T for drought studies in Gumbel’s technique 
is written as: 


Qt 


Q 


i + 


c v {y T ~ y n ) 


( 9 ) 


Considering a stream flow drought index (£D/ m ) with a rank 
m that corresponds to a particular stream flow Q then by 
applying Equation (9) this results to: 


Generally the results show that, the minimum and maximum 
drought episodes were exhibited in gauge stations 4AB05 
(Amboni) and 4CC03 (Yatta furrow) with absolute SDI rang¬ 
ing from 0.667 tol.265 and 1.213 to 2.42 for2 and 1000-year 
return period respectively. 


Q m = Q 


i + 



( 10 ) 


When Equations (9) and (10) are reorganized and then divid¬ 
ing Equation (9) by (10) the relation becomes: 



“♦— T5 — »-T10 —A— T20 -X-T50 ‘ T100 

Figure 4: Q T lQ versus C v for 4ED01 (Kamburu) gauge sta¬ 
tion in upper Tana River basin 


Qm-Q y,n-y n (11) 

Simplifying Equation (11) leads to: 

Qt ~Q _ TV ~ y„ (i 2 ) 

Q m -Q y m - y n 

From Gumbel’s method the estimated values of y n as given 
in Table 5 are 0.5236 and 0.5745 for the data record of 20 
and 1000 years respectively. Thus the parameter y n can be 
assumed to be a constant that arbitrarily lies between 0.5236 
and 0.5747 and estimated using the relation: 
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Kamburu 4ED01 1.112 19.463 1.319 16.182 1.453 14.071 1.592 12.044 1.764 9.332 1.893 7.382 2.022 5.401 2.192 2.447 
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y„ =0.5236 + ^x (0.5747 -0.5236) = 0.54915 = 0.5 (13 ) 

The computed value of 0.55 corresponds to hydrological 
drought event of a particular severity with a return period of 
50 years (Table 5). 


Table 5: Typical values of y n and o n in Gumbel’s 
method 



20 

0.5236 

1.0628 

70 

0.5548 

1.1824 

25 

0.5309 

1.0915 

75 

0.5559 

1.1898 

30 

0.5362 

1.1124 

80 

0.5569 

1.1938 

35 

0.5403 

1.1285 

85 

0.5578 

1.1973 

40 

0.5436 

1.1413 

90 

0.5589 

1.2007 

45 

0.5463 

1.1518 

95 

0.5593 

1.2038 

50 

0.5465 

1.1607 

100 

0.5600 

1.2065 

55 

0.5504 

1.1681 

200 

0.5672 

1.2359 

60 

0.5521 

1.1747 

500 

0.5724 

1.2588 

65 

0.5536 

1.1808 

1000 

0.5745 

1.2685 


N = number of years, y n =expected mean of the data, 
o n =expected standard deviation of the data 

Source: Gumbel (1958) 


Based on Gumbel’s method the return period T and 
y T , are determined using the following functions: 


T = 


N + l 


(14) 


y r =-lnln^y (15) 

For a particular stream flow drought index corresponding to 
stream flow O with a rank in, the value of y is detennined 
from: 



This relation reduces to: 


y,n = - In In 


' N + l ) 
K N +1 — m J 


(16) 


(17) 


Qt-Q 

- In In | 

T 1 0.55 
yr-ij 

^ m= 1 

42 >H 

J m= 1 

f ^ + 1 1 055 

^ N +1 — m ) 


( 19 ) 


From the principles of Schulz (1973), it can be shown that: 

- In In f —1 = In T — f ——I-j- + —1 (20) 

(r-l) {2 T 24 T 2 ST 3 ) v ' 

Neglecting the quantities inside the brackets on the right side 
of Equation (20), it reduces to: 


-lnln 


< J > 


= ln T 


( 21 ) 


Also neglecting the quantities in the brackets for Equation 
(17) yields: 


y m = - In In 


N + l 


A , fAr+n 

= In - 

V m 


( 22 ) 


^N + l — rri; 

Therefore, Equation (19) can be reduced further into the fol¬ 
lowing relation: 


Qt~Q 

\YjQm~Q 

4 m=l 


1 

3 


In T- 0.55 


X ln 


N + l 


0.55 


(23) 


Table 4 show stream flow values corresponding to absolute 
stream flow drought index that were computed for different 
return periods using Equations (22) and (23) based on data 
acquired for the upper Tana River basin. 


By redefining Equation (23) in X and Y as shown in Equa¬ 
tions (24) and (25), and plotting the corresponding data for 
the upper Tana River basin gives Figure 5. 


X = 


lnT-0.55 


1, I 

'N + l') 

— In 

0.55 

3 

, m ) 


(24) 


Y= Q-Qt 

e-)ie„ 

^ m =1 


(25) 


Substituting Equation (15) into Equation (10) results to: 


Qt~Q 

Q m -Q 



(18) 


The results from Figure 5 show that all the gauged stations 
exhibit strong linear correlation. 
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a) x b) 


Figure 5: Fitted Data for 4BC02 (Tana Sagana) and 4AC03 
(Sagana) gauged stations 

The plot shows that there is a strong correlation between 
the expressions on the left and right side of Equation (23), 
with correlation coefficients at gauge stations of IDs 4BC02 
(Tana Sagana) and 4AC03 (Sagana) of 0.826 and 0.793 re¬ 
spectively. This means that the stream flow Q of any return 
period can be determined from mean flows of the annual 
minimum and the average of the first three minimum mean 
stream flows for the upper Tana river basin as per Equation 
(23). This is found to be consistence with the similar plot de¬ 
veloped for flood estimation by Al-Mashindani et al. (1978). 


CONCLUSION 

i) From the study SDI, stream flows have been explored 
and their corresponding return periods estimated. It is 
concluded that the computed absolute SDI vary across 
the gauge stations and increase while the correspond¬ 
ing stream flow decline with increase in return period 

ii) A simplified mathematical model for estimating hy¬ 
drological drought event that uses the mean of the an¬ 
nual minimum and average of the first three minimum 
stream flows as input variables is developed for differ¬ 
ent return periods in the upper Tana River basin. 

iii) Critical points upon which design of water storage 
systems can be based are identified for different gauge 
stations. These indicate the level of hydrological 
drought beyond which the water facilities is signifi¬ 
cantly affected by the drought. 
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